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Abstract- The  fourth  passive  circuit  element,  memristor,  has 
attracted  increased  attentions  since  the  first  real  device  was 
discovered  by  HP  Lab  in  2008.  Its  distinctive  characteristic  to 
record  the  historic  profile  of  the  voltage/current  through  itself 
creates  great  potentials  in  future  system  design.  However,  as  a 
nano-scale  device,  memristor  is  facing  great  challenge  on 
process  variation  control  in  the  manufacturing.  In  this  work,  we 
analyze  the  impact  of  the  geometry  variations  on  the  electrical 
properties  of  both  Ti02  thin-film  and  spintronic  memristors, 
including  line  edge  roughness  and  thickness  fluctuation.  A 
simple  algorithm  was  proposed  to  generate  a  large  volume  of 
geometry  variation-aware  three-dimensional  device  structures 
for  Monte-Carlo  simulations.  Our  simulation  results  show  that 
due  to  the  different  physical  mechanisms,  Ti02  thin-film 
memristor  and  spintronic  memristor  demonstrate  very  different 
electrical  characteristics  even  when  exposing  them  to  the  same 
excitations  and  under  the  same  process  variation  conditions. 

I.  Introduction 

Nearly  forty  years  ago,  Professor  Chua  predicted  the 
existence  of  memristor  -  the  fourth  fundamental  circuit 
element,  to  complete  the  set  of  passive  devices  that 
previously  includes  only  resistor,  capacitor,  and  inductor  [1]. 
In  2008,  the  first  physical  realization  of  memristor  was 
demonstrated  by  HP  Lab,  in  which  the  memristive  effect  was 
achieved  by  moving  the  doping  front  along  Ti02  thin-film 
device  [2].  Soon,  memristive  systems  on  spintronic  devices 
were  proposed  [3]. 

The  unique  properties  of  memristor  create  great 
opportunities  in  future  system  design.  For  instance,  the  non¬ 
volatility  and  excellent  scalability  make  it  a  promising 
candidate  as  the  next-generation  high-performance  high- 
density  storage  technology  [4].  The  applications  of 
memristive  behavior  [5]  in  electronic  neural  network  also 
have  been  extensively  studied  [6]  [7]. 

As  process  technology  shrinks  down  to  decananometer 
(sub-5 Onm)  scale,  device  parameter  fluctuations  incurred  by 
process  variations  have  become  a  critical  issue  affecting  the 
electrical  characteristics  of  devices  [8].  The  situation  in  a 
memristive  system  could  be  even  worse  when  utilizing  the 
analog  states  of  the  memristor  in  design.  Previous  works  on 
memristor  variation  analysis  mainly  focused  on  its  impacts 
on  non-volatile  memory  design  [4]  [9].  However,  the 


systematic  analysis  and  quantitative  evaluation  on  how 
process  variations  affect  the  memristive  behaviors  are  still 
missing.  Our  work  explores  the  implications  of  the  device 
parameters  of  memristors  to  the  circuit  design  by  taking  into 
account  the  impacts  of  process  variations.  The  evaluations 
were  conducted  based  on  both  theoretical  analysis  and  Monte 
Carlo  simulations. 

The  device  geometry  variations  significantly  influence  the 
electrical  properties  of  nano-devices  [10].  For  example,  the 
random  uncertainties  in  lithography  and  patterning  processes 
lead  to  the  random  deviations  of  line  edge  print-images  from 
its  ideal  pattern,  which  is  called  line  edge  roughness  (LER) 
[11].  Thickness  fluctuation  is  caused  by  deposition  process  in 
which  mounds  of  atoms  form  and  coarsen  over  time.  As 
technology  shrinks,  the  geometry  variations  do  not  decrease 
accordingly.  In  this  work,  we  propose  an  algorithm  to 
generate  a  large  volume  of  three-dimensional  (3D)  memristor 
structures  to  mimic  the  geometry  variations.  The  LER  model 
is  based  on  the  latest  LER  characterization  method  for 
electron  beam  lithography  (EBL)  technology  from  top-down 
scanning  electron  microscope  (SEM)  measurement  [12]. 
Other  process  variations  such  as  random  discrete  doping 
(RDD)  could  also  result  in  the  fluctuations  of  the  electrical 
properties  of  devices.  However,  because  the  existing 
memristors  are  all  based  on  the  thin  film  deposition 
technology,  the  local  randomness  of  RDD  is  not  as 
significant  as  geometry  variations,  and  hence,  is  not  covered 
in  this  work. 

Memristive  function  can  be  achieved  by  various  materials 
and  device  structures.  The  impact  of  the  process  variations  on 
the  electrical  properties  of  different  memristors  could  be  very 
different  even  under  the  same  excitations.  Two  types  of 
devices,  Ti02  thin-film  memristor  [3]  and  spintronic 
memristor  [13],  are  analyzed  and  evaluated  in  this  work, 
respectively. 

The  organization  of  this  paper  is  as  follows:  Section  II 
introduces  the  fundamentals  of  memristor  theory  and  the 
physical  mechanisms  of  Ti02  thin-film  memristor  and 
spintronic  memristor;  Section  III  analyzes  the  memristor 
models  under  process  variations;  Section  IV  explains  our 
proposed  algorithm  for  3D  memristor  structure  generation; 
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Figure  1.  Ti02  thin-film  memristor.  (a)  structure,  and  (b)  equivalent 
circuit. 


Figure  2.  TMR-based  spintronic  memristor.  (a)  structure,  and  (b) 
equivalent  circuit. 


Section  V  presents  and  analyzes  the  simulation  results  of  two 
types  of  memristors;  and  Section  VI  concludes  our  work. 

II.  Preliminaries 


A.  Memristor  Theory 

The  original  definition  of  memristor  is  derived  from  the 
completeness  of  circuit  theory:  besides  resistor,  capacitor  and 
inductor,  there  must  exist  the  fourth  basic  two-terminal 
element  that  uniquely  defines  the  relationship  between  the 
magnetic  flux  ( cp )  and  the  electric  charge  ( q )  passing  through 
the  device  [1],  or 

dcp  =  M  •  dq.  (1) 

Considering  that  magnetic  flux  and  electric  charge  are  the 
integrals  of  voltage  and  current  over  time,  respectively,  the 
definition  of  memristor  can  be  generalized  as: 

(P  =  Af  (*>,/)•/ 

{da)/dt  =  { 

Here,  m  is  a  state  variable;  M(m,  /)  represents  the 
instantaneous  memristance,  which  varies  over  time.  For  a 
“pure”  memristor,  neither  M(m,  /)  nor  is  an  explicit 

function  of  /  [5]. 


B.  T1O2  Thin-film  Memristor 


In  2008,  HP  Lab  demonstrated  the  first  memristive  device 
in  Ti02  thin-film  [2].  Figure  1  illustrates  the  conceptual  view 
and  the  corresponding  variable  resistor  model,  which  is 
equivalent  to  two  series-connected  resistors.  Here,  RL  (Rfi  is 
used  to  denote  the  lowest  (highest)  resistance.  The  overall 
resistance  can  be  expressed  as 

M(a)  =  a-RL  +  (l-d)-RH.  (3) 


a  (0  <  a  <  1)  is  the  relative  doping  front  position, 
which  is  the  ratio  of  doping  front  position  over  the  total 
thickness  of  Ti02  thin-film.  The  velocity  of  doping  front 
movement  v(t),  which  is  driven  by  the  voltage  applied  across 
the  memristor  P(t)  can  be  expressed  as 


Rl 


no 

M(a)’ 


(4) 


where,  fiv  is  the  equivalent  mobility  of  dopants;  h  is  the 
total  thickness  of  the  Ti02  thin  film;  and  M(a)  is  the  total 
memristance  when  the  relative  doping  front  position  is  a. 


structure  [3] [13].  We  choose  tunneling  magneto-resistance 
(TMR)-based  structure  illustrated  in  Figure  2(a)  as  the 
objective  of  this  work.  The  lowest  (highest)  resistance  is 
denoted  as  RL  (Rn)- 

As  shown  in  Figure  2(b),  the  overall  resistance  of  a  TMR- 
base  spintronic  memristor  can  be  modeled  as  two  parallel- 
connected  resistors  with  resistances  RL/a  andPH/(l  — a ), 
respectively  [13].  Here  a  (0  <  a  <  1)  represents  the  relative 
domain  wall  position  as  the  ratio  of  the  domain  wall 
position  (x)  over  the  total  length  of  the  free  layer  (L).  The 
overall  memristance  can  be  expressed  as 


Mid)  = 


Rh'Rl 

RH-a+RL-(l-a)' 


(5) 


The  domain-wall  velocity  v(t)  is  proportional  to  the 
current  density  /  [14].  We  have 


m  = 


no 

M(a)-L-z 


,  and 


(6) 


/<£■  <7) 

Here  Tv  is  the  domain  wall  velocity  coefficient,  which  is 
related  to  device  structure  and  material  property.  L  and  z  are 
the  total  length  and  width  of  the  spintronic  memristor, 
respectively.  The  domain  wall  movement  in  a  spintronic 
memristor  happens  only  when  the  current  density  ]{t)  is 
above  the  critical  current  density  (Jcr )  [14]  [15]. 


III.  Statistical  Analysis 

The  actual  length  ( L )  and  width  (z)  of  a  memristor  is 
affected  by  line  edge  roughness  (LER).  The  variation  of 
thickness  ( h )  of  a  thin  film  structure  is  usually  described  by 
thickness  fluctuations.  As  a  matter  of  convenience,  we  define 
that,  the  impact  of  process  variations  on  any  given  variable 
can  be  expressed  as  a  factor  9  =  cn'/co,  where  a)  is  the  ideal 
value,  and  a)'  is  the  actual  value  under  process  variations. 
The  ideal  geometry  dimensions  of  Ti02  thin-film  memristor 
and  spintronic  memristor  used  in  this  work  are  summarized 
in  Table  I  [2]  [13]. 

A.  Ti02  Thin-Film  Memristor 

In  a  Ti02  thin-film  memristor,  the  current  passes  through 


C.  Spintronic  Memristor 

Among  all  the  spintronic  memristive  devices  proposed 
recently,  the  one  based  on  magnetic  tunneling  junction  (MTJ) 
could  be  the  most  promising  one  because  of  its  simple 


Table  I.  The  Device  Dimensions  of  Memristors. 


Length  ( L ) 

Width  (z) 

Thickness  (h) 

Thin-film 

50  nm 

50  nm 

10  nm 

Spintronic 

200  nm 

10  nm 

7  nm 
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the  device  along  thickness  ( h )  direction.  Ideally  the  doping 
front  has  an  area  S=L  z.  To  simulate  the  impact  of  LER  on 
the  electrical  properties,  the  memristor  device  is  divided  into 
many  small  filaments  between  the  two  electrodes.  Each 
filament  z  has  a  cross-section  area  ds  and  a  thickness  h.  The 
ideal  upper  bound  and  lower  bound  of  the  memristance  for 
filament  z  are 


RiiH=RH-j-,mdRt,L  =  RL-j-. 


(8) 


For  filament  z,  the  actual  upper  and  the  lower  bound 
under  the  process  variations  can  be  expressed  as 

Kh  =  Ri,H  ■  and  R[l  =  Rul  •  ft.  (9) 

Here  #1?s  and  6^  represent  the  variations  caused  by  2-D 
LER  and  thickness  fluctuations,  respectively. 

The  doping  front  velocity  in  filament  z  by  considering 
process  variations  can  be  calculated  as: 

rI,l  no 


(10) 


m'iWY 

Here  hr  and  M[  are  the  actual  thickness  and  memristance 
of  filament  z.  Then,  we  can  get  a  set  of  related  equations  for 
filament  z,  including  the  doping  front  position 
rt. 


a'i(t)  =  f0  v[(t)  ■  dr, 

the  corresponding  memristance 

Ml(a'i)  =  a'rR'u  +  (l-a'i)-R'i'H, 
and  the  current  through  the  filament  z 

/'(t)  = 


(11) 


(12) 


(13) 


By  combining  Eq.  (10)  -  (13),  the  doping  front  position 
in  every  filament  z  under  process  variations  a[{t)  can  be 
obtained  by  solving  the  differential  equation 


da'jjt)  _  R[l _ V(t) _ 

dt  ~{iv'  h'2  ' 


(14) 


Eq.  (14)  indicates  that  the  behavior  of  the  doping  front 
movement  is  dependent  on  the  specific  electrical  excitations, 
i.e.,  V(t).  For  instance,  if  applying  a  sinusoidal  voltage  source 
such  as 

V  =  Vm  •  sin  (2nf  •  t),  (15) 

the  corresponding  doping  front  position  of  filament  z  can  be 
expressed  as: 


a[(t ) 


RfH+A-B(t)-£—2C[l]-A~-^ 
ei.h  0i.h 


Where,  A  =  RiH  -  RiL ,  B(t )  =  rv  •  Rt 


0  <  a[  <  1.  (16) 

COS(27T/-t) 


f'h2'U 


and  C[l]  is  an  initial  state  constant. 

The  existence  of  term  B(t)  stands  for  the  impact  of 
electrical  excitation  on  doping  front  position.  The  terms  0{  s 
and  h  stand  for  the  impact  of  both  LER  and  thickness 
fluctuations  on  memristive  behavior,  respectively.  Moreover, 
the  impact  of  the  geometry  variations  on  the  electrical 
properties  of  memristors  could  be  affected  by  the  electrical 
excitations.  For  example,  we  can  set  a( 0)  =  0  to  represent 
the  case  that  the  Ti02  memristor  starts  from  M(0)=R H.  In 


such  a  condition,  the  doping  front  position  a[{t)  can  be 
calculated  as: 


Ri,H~ 


am 


9i,h  ei,h 


0  <  a[  <  1, 


(17) 


which  is  affected  only  by  thickness  fluctuations  and  electrical 
excitations.  LER  won’t  disturb  a[{t)  if  the  Ti02  thin-film 
memristor  has  an  initial  state  a(0)  =  0. 

The  overall  memristance  of  the  memristor  can  be 
calculated  as  the  total  resistance  of  all  n  filaments  connected 
in  parallel.  When  n  increases  to  oo?  we  have 

1  ,  r-i !  i 


R'h  =  R* 


XH  '  Cei*/e‘,s-dl’  and  Ri  Ri  '  Ce‘*/Btsdl  ' 

The  overall  current  through  the  memristor  is  the  sum  of 


(18) 


the  current  through  each  filament: 


(19) 


The  instantaneous  memristance  of  the  overall  memristor 
can  be  defined  as 


M'(t) 


no 

'  i'it ) 


(20) 


CVMl-di 

Considering  that  the  doping  front  movement  in  each 
filament  won’t  be  the  same  because  h[  varies  due  to  thickness 
fluctuation  (and/or  the  roughness  of  the  electrode  contact), 
we  define  the  average  doping  front  position  of  the  whole 
memristor  as: 

R'H-M'(t ) 


ce'(t)  = : 


R'u-R'j 


(21) 


B.  Spintronic  Memristor 

Because  the  length  of  a  spintronic  memristor  is  usually 
much  longer  than  the  other  dimensions,  the  impact  of  the 
variance  in  length  on  the  device’s  electrical  properties  can  be 
ignored.  In  our  analysis,  the  device  can  be  chopped  into 
infinite  segments  in  the  length  direction.  For  a  segment  z,  the 
upper  and  lower  bounds  of  memristance  are 


R'i,H  =  ri,h  R'liL  =  Rl,L-e-*. 


(22) 


, - 

“ l,Z  “lz 

Here  we  assume  the  ideal  memristance  changes  linearly 
within  the  domain  wall,  or  R ijH  changes  linearly  from  7?jjL  to 
^k,H  when  j  <  i  <  k.  Here  j  and  k  are  the  two  segments  at  the 
two  boundaries  of  domain  wall  and  connected  to  the 
magnetic  domains  with  either  the  low  or  the  high  resistance 
states.  The  memristance  of  each  segment  is 

(RU'i<<*' 

Wi,H>  i  >  «r 

So  for  overall  resistance  R'H  and  R'h,  we  have 

1  1  n  /  n  1 


Mj(a')  = 


(23) 


R'„-- 


A.JJ  '  roo  ,  and  Rj  —  R j  m  roo 

H  J0  0i,z/0i,h-di  1  J0  eiiZ/eUh-di 

Then  the  memristance  of  the  whole  device  is 
1  1 


"'(«')  =  7^T 


tf-j-dCO+tij-dCO  f0° 


Ri,L°i,h 


1  Ri,Hdi,h 


d{  i) 


(24) 


-•  (25) 


Here  the  width  of  every  segments  r,  varies  segment  by 
segment  due  to  the  LER  effect.  The  statistical  behavior  of 
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Figure  4.  An  example  of  3-D  memristor  structure. 


Figure  3.  The  flow  of  3D  memristor  structure  generation  including 
geometry  variations. 

spintronic  memristors  can  still  be  evaluated  by  Monte-Carlo 
simulation  in  Section  V. 

We  assume  the  current  density  applied  on  the  domain 
wall  J'(t)  is  the  one  of  the  segments  i  in  which  the  middle  of 
domain  wall  located: 


nt)=n 


V(t) 


(26) 


Then  the  domain  wall  velocity  under  process  variations 
can  be  calculated  as: 


=  i.„  =  f0; 


\r,  r>Ja 

r  <Ja 


(27) 


IV.  3D  Memristor  Structure  Modeling 

Analytic  modeling  is  a  fast  way  to  estimate  the  impacts  of 
process  variations  on  memristors.  However,  we  noticed  that 
modeling  some  variations  analytically  e.g.  LER  may  be 
beyond  the  capability  of  analytic  model  [12].  The  silicon  data 
of  these  variations,  however,  is  usually  very  hard  to  obtain 
simply  due  to  intelligent  property  protection.  To  improve  the 
accuracy  of  our  evaluations,  we  create  a  simulation  flow  to 
generate  3D  memristor  samples  with  the  geometry  variations 
including  LER  and  thickness  fluctuation.  The  correlation 
between  the  generated  samples  and  the  real  silicon  data  are 
guaranteed  by  the  sanity  check  of  the  LER  characterization 
parameters.  The  flow  is  shown  in  Figure  3. 

Many  factors  affecting  the  quality  of  the  line  edges  show 
different  random  effects.  Usually  statistical  parameters  such 
as  auto-correlation  function  (ACF)  and  power  spectral 


Table  II.  The  Parameters/Constraints  in  LER 


Characterization. 


Parameters 

Constraints 

^LF 

0.8  nm 

<TlER 

2.5nm  -  3.5nm 

/nax 

1.8  MHz 

<TlWR 

4.0nm~  5.0  nm 

Euf 

0.4  nm 

skewness 

0.1  ~0.2 

/ 

/ 

kurtosis 

2.9-3. 1 

Table  III.  Memristor  Device  and  Electrical  Parameters. 


Ti02  thin-fiim  memristor  [2 

Rl  m 

Ru  (-O) 

My  (mW'1) 

/ 

rn(V) 

/(Hz) 

100 

16000 

/ 

1 

0.5 

Spintronic  memristor  [13] 

Rl  m 

Rk(Q.) 

rv(nm3-C-') 

Jcr  (A/nm2) 

vm  (V) 

/(Hz) 

2500 

7500 

2.01xl016 

2.00x1 0'8 

2 

10M 

density  (PSD)  are  used  to  describe  the  property  of  the  line 
edges.  Considering  that  LER  issue  is  related  to  fabrication 
process,  we  mainly  target  the  nano-scale  pattern  fabricated  by 
electron  beam  lithography  (EBL).  The  measurements  show 
that  under  such  a  condition,  the  line  edge  profile  has  two 
important  properties:  (1)  the  line  edge  profile  demonstrates 
regular  oscillations,  which  are  caused  by  periodic 
composition  in  EBL  fabrication  system;  and  (2)  the  line  edge 
roughness  mainly  concentrates  in  a  low  frequency  zone  [12]. 

To  generate  line  edge  samples  close  to  the  real  cases,  we 
can  equally  divide  the  whole  line  edge  into  many  segments, 
say,  n  segments.  The  LER  of  the  ith  segment  can  be  modeled 
by 

LERi  =  Lip  •  sin (fmax  •  x{)  +  LHF  •  p^.  (28) 

The  first  term  in  the  right  side  of  Eq.  (28)  represents  the 
regular  disturbance  at  the  low  frequency  range,  which  is 
modeled  as  a  sinusoid  function  with  amplitude  CLF./max  is  the 
mean  of  the  low  frequency  range  derived  from  PSD  analysis. 
Without  loss  of  generality,  an  uniform  distribution  Xi  <e  U(- 
1,1)  is  used  to  represent  an  equal  distribution  of  all  frequency 
components  in  the  low  frequency  range.  The  disturbance 
from  high  frequency  domain  is  also  taken  into  account  by  the 
second  term  on  the  right  side  of  Eq.  (28)  as  a  Gaussian  white 
noise  with  amplitude  f-  Here  p{  follows  the  normal 
distribution  N (0,1)  [12].  The  actual  values  of  ZLF,  ZH f  and 
/max  are  determined  by  ACF  and  PSD. 

To  ensure  the  correlation  between  the  generated  line  edge 
samples  with  the  measurement  results,  we  introduce  four 
constraints  to  conduct  sanity  check  of  the  generated  samples: 

•  0ler-  the  root  mean  square  (RMS)  of  LER; 

•  ^lwr:  the  RMS  of  line  width  roughness  (LWR); 

•  skewness:  the  symmetry  of  the  amplitude  of  the  line 
edge;  and 

•  kurtosis:  the  steepness  of  the  amplitude  distribution 
curve. 

The  above  four  parameters  are  widely  received  in  LER 
characterization  and  can  be  achieved  from  measurement 
results  directly.  Only  the  line  edge  samples  that  satisfy  the 
constraints  will  be  taken  for  the  Monte-Carlo  simulations. 
Table  II  summarizes  the  parameters  used  in  our  algorithm, 
which  are  correlated  with  the  characterization  method  and 
experimental  results  in  [12]. 

The  thickness  fluctuation  is  caused  by  the  random 
uncertainties  in  sputter  deposition  or  atomic  layer  deposition. 
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Table  IV.  3  a  min./max.  of  Ti02  Memristor  Parameters. 


Sinusoidal 

Voltage 

LER  only 

Thickness  only 

Overall 

-3g 

+3  a 

-3o 

+3  a 

-3  a 

+3  a 

-^H  & 

-5.4% 

4.1% 

-5.5% 

4.8% 

-6.4% 

7.3% 

M(a) 

-5.4% 

4.1% 

-37.1% 

20.8% 

-36.5% 

24.1% 

aft) 

0.0% 

0.0% 

-13.3% 

27.5% 

-14.7% 

27.4% 

v(a) 

0.0% 

0.0% 

-9.3% 

15.6% 

-10.4% 

16.9% 

i(a) 

-4.7% 

5.7% 

-9.3% 

15.7% 

-10.7% 

17.2% 

Power 

-4.7% 

5.7% 

-8.8% 

14.1% 

-10.1% 

15.6% 

It  has  a  relatively  smaller  impact  than  LER  and  can  be 
modeled  as  Gaussian  distribution.  Considering  that  the 
memristors  in  this  work  have  relatively  bigger  dimensions  in 
horizontal  plane  than  thickness  direction  (shown  in  Table  I), 
we  also  considered  roughness  of  electrode  contact  in  our 
simulation:  The  means  of  the  thickness  of  every  memristors 
are  generated  by  assuming  it  follows  the  Gaussian 
distribution  among  all  simulated  memristors.  Then  every 
memristor  is  divided  into  many  filaments  between  the  two 
electrodes.  The  roughness  of  electrode  contracts  is  modeled 
as  the  disturbance  on  the  thickness  of  each  filament.  We 
assume  that  both  thickness  fluctuations  and  electrode  contact 
roughness  follow  Gaussian  distributions  with  a  deviation 
g=2%  of  thin  film  thickness. 

Figure  4  is  an  example  of  3-D  spintronic  memristor 
structure  generated  by  our  flow. 

V.  Experimental  Results 

A.  Simulation  Setup 

To  evaluate  the  impacts  of  process  variations  on  the 
electrical  properties  of  memristors,  we  conducted  Monte- 
Carlo  simulations  with  10,000  qualified  3-D  device  samples 
generated  by  our  proposed  flow.  A  sinusoidal  voltage  source 
in  Eq.  (15)  is  applied  as  excitation.  The  initial  state  of 
memristor  is  set  as  M(cf=0)=Ry{.  The  device  and  electrical 
parameters  used  in  our  simulations  are  summarized  in  TABLE 
III.  Both  separated  and  combined  effects  of  geometry 
variations  on  various  electrical  properties  are  analyzed, 
including: 

•  the  distribution  of  and  RL ; 

•  the  change  of  memristance  M(t)  and  M(a); 

•  the  velocity  of  wall  movement  v(a); 

•  the  current  through  memristor  i(t );  and 

•  the  I-V  characteristics. 

B.  HO2  Thin-Film  Memristor 

For  Ti02  thin-film  memristors,  the  ±3 a  (minimal/ 
maximal)  values  of  the  device/electrical  parameters  in  the 
percentage  of  the  corresponding  ideal  values  are  summarized 
in  Table  IV.  For  those  parameters  that  vary  over  time,  we 
consider  the  variation  at  each  time  step  of  all  the  devices.  The 
simulation  results  considering  only  either  FER  or  thickness 
fluctuations  are  also  listed.  To  visually  demonstrate  the 
overall  impact  of  process  variations  on  the  memristive 
behavior,  the  dynamic  responses  of  100  Monte  Carlo 
simulations  are  shown  in  Figure  5. 


Figure  5.  Simulation  results  for  Ti02  thin-film  memristors.  The  blue 
curves  are  from  100  Monte-Carlo  simulations,  and  red  lines  are  the 
ideal  condition.  From  top  left  to  right  bottom,  the  figures  are  Ru  vs. 
Rl;  M(t)  vs.  /;  v  vs.  a ;  a  vs.  /;  /  vs.  /;  and  I-V  characteristics. 

Table  IV  shows  that  the  static  behavior  parameters,  i.e., 
Rn  and  R^,  are  affected  in  a  similar  range  by  FER  and 
thickness  fluctuations.  This  is  consistent  to  the  analytical 
results  in  Eq.  (18):  0S  and  6^  have  the  similar  weights  on  the 
variation  of  R'H  and  R[ . 

However,  thickness  fluctuations  show  a  much  more 
significant  impact  on  the  memristive  behaviors  such  as  v(t), 
a(t)  and  M(a),  than  FER  does  because  the  doping  front 
movement  is  along  the  thickness  direction:  vft)  is  inversely 
proportional  to  the  square  of  thickness,  and  aft)  is  the  integral 
of  vft)  over  time  as  shown  in  Eq.  (10)  and  (11).  Due  to  the 
similar  reason,  thickness  fluctuations  significantly  affect  the 
instantaneous  memristance  M(a)  as  well. 

Since  the  thickness  of  Ti02  thin-film  memristor  is 
relatively  small  compared  to  other  dimensions,  we  assume 
the  doping  front  cross-section  area  is  a  constant  along 
thickness  direction  in  our  simulation.  The  impacts  of  FER  on 
aft)  or  vft)  are  negligible  compared  to  that  of  the  thickness 
fluctuation  as  shown  in  TABLE  IV. 

An  interesting  observation  in  Figure  5  is  that  as  the 
doping  front  a  moves  toward  1,  the  velocity  v  regularly 
grows  larger  and  reaches  its  peak  at  the  half  period  of  the 
sinusoidal  excitation,  i.e.  t=  Is.  This  can  be  explained  by  Eq. 
(16):  the  memristance  is  getting  smaller  as  a  moves  toward  1. 
With  the  same  input  amplitude,  a  smaller  resistance  will 
result  in  a  bigger  current  and  therefore  a  bigger  variation  on 
vft).  Similarly,  memristance  M(a)  reaches  its  peak  variance 
when  a  is  close  to  1 . 

C.  Spintronic  Memristor 
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Table  V.  3  a  min./max.  of  Spintronic  Memristor  Parameters. 


Sinusoidal 

Voltage 

LER  only 

Thickness  only 

Overall 

-3  a 

+3  a 

-3  a 

+3  a 

-3  a 

+3  a 

-15.3% 

22.9% 

-6.1% 

5.8% 

-16.4% 

20.9% 

M(a) 

-15.1% 

23.3% 

-11.0% 

11.0% 

-16.3% 

21.1% 

aft) 

-9.7% 

8.1% 

-8.4% 

9.5% 

-11.8% 

8.1% 

v(a) 

-10.7% 

22.1% 

-9.1% 

9.9% 

-21.5% 

22.5% 

i(a) 

-18.5% 

18.5% 

-8.9% 

10.1% 

-17.7% 

17.8% 

Power 

-18.4% 

18.6% 

-8.3% 

9.4% 

-17.8% 

17.8% 

The  ±3 a  values  of  the  device/electrical  parameters  based 
on  10,000  Monte-Carlo  simulations  for  spintronic  memristors 
are  summarized  in  TABLE  V.  The  visual  demonstration  of  100 
Monte-Carlo  simulations  with  a  sinusoidal  voltage  excitation 
is  shown  in  Figure  6. 

In  spintronic  memristors,  the  impact  of  LER  is  more  than 
that  of  thickness  fluctuations.  This  is  because  the  direction  of 
the  domain  wall  movement  is  perpendicular  to  the  direction 
of  spin-polarized  current.  The  impacts  of  thickness 
fluctuations  on  very  small  segments  cancel  each  other  during 
the  integral  along  the  direction  of  the  domain  wall 
movement. 

“LER  only”  simulation  results  show  more  impact  on  the 
electrical  properties  at  +3 a  corner  than  that  at  -3a  corner. 
This  is  because  the  variation  of  line  width  is  the  dominant 
factor  of  the  variation  of  electrical  properties  of  spintronic 
memristors  and  the  line  edge  profiles  used  in  our  LER 
parameters  have  a  right-biased  feature  [12].  Since  normal 
distribution  is  assumed  for  the  variations  of  thickness,  6^  has 
approximately  symmetric  impacts  on  ±3 a  corners 

The  impact  of  LER  on  the  memristive  parameters  v(t), 
aft)  and  M(a)  is  also  larger  than  thickness  variations.  Again, 
the  impacts  of  thickness  fluctuations  on  very  small  segments 
cancel  each  other  during  the  integral  along  the  direction  of 
the  domain  wall  movement. 

VI.  Conclusion 

In  this  work,  we  evaluate  the  impact  of  different  geometry 
variations  on  the  electrical  properties  of  two  different  types 
of  memristors  -  Ti02  thin-film  and  spintronic  memristors,  by 
conducting  the  analytic  modeling  analysis  and  Monte-Carlo 
simulations.  We  investigate  the  different  responses  of  the 
static  and  memristive  parameters  under  various  process 
variations  and  explore  its  implication  to  the  electrical 
properties  of  the  memristors.  A  simple  LER  sample 
generation  algorithm  is  proposed  to  speed  up  the  related 
Monte-Carlo  simulations.  To  author’s  best  knowledge,  this  is 
the  first  work  that  conduct  the  quantitative  evaluation  and 
comparison  on  the  impact  of  geometry  variations  to  the  two 
types  of  memristors. 
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